Non-uniform thermal magnetization noise in thin films: application to GMR heads 
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A general scheme is developed to analyze the effect of non- 
uniform thermal magnetization fluctuations in a thin film. 
The normal mode formalism is utilized to calculate random 
magnetization fluctuations. The magnetization noise is pro- 
portional to the temperature and inversely proportional to the 
film volume. The total noise power is the sum of normal mode 
spectral noises and mainly determined by spin-wave standing 
modes with an odd number of oscillations. The effect rapidly 
decreases with increasing mode number. An exact analytical 
calcutaion is presented for a two-cell model. 



Our approach is purely classical and final result does not 
contain fi. 

The energy of the spin system can be written as 



where 
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dd 
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I. INTRODUCTION 



is the energy of anisotropy, is the unit vector along 
the axis of in-plane uniaxial anisotropy, H^^- is the lo- 
cal uniaxial anisotropy field, H^^^ is the local hard-axis 
anisotropy field. 
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Recently thermal magnetization fluctuations have been 
shown to be an important issue for thin film magnetore- 
sistive heads used for high-density magnetic recording [1] . 
Noise estimates have been initially made for the simplest 
model of a uniformly magnetized film [1-3]. In order to 
develop a more realistic model it is necessary to analyze 
the role of non- uniform magnetization fluctuations. In 
this paper we propose a general scheme to calculate the 
effect of non-uniform thermal magnetization fluctuations 
in an ultrathin film. After discretizing the film into mag- 
netic cells, the well-known normal mode formalism is uti- 
lized. By doing so we represent the magnetization noise 
in terms of independent standing spin waves subjected to 
white thermal noise. As an example, the explicit solution 
is given for a two-mode model. 



II. CALCULATION SCHEME 



Let us consider a thin film of the volume V 



N 



dr, where 



< dy, dz 



are the thickness, width 



and length, respectively. The thickness is assumed to be 
of the order or less than exchange length {d^ < VA/Ms) 
and therefore non-uniform magnetization dynamics along 
the X direction will be neglected. The film is discrctized 
in the yz plane into N identical magnetic cells. Each cell 
is assumed to be uniformly magnetized and can be rep- 
resented as a classical spin Sj — —'M.{rj)Vc/hj, where 
M(rj) is the magnetization in the cell (|M(rj)| = Mg), 
Vc = V/N is the cell volume, Ms is the saturation mag- 
netization and 7 is the gyromagnetic ratio. Planck's con- 
stant h is used for convenience as a dimensional variable. 



(3) 



is the Zeeman energy with local magnetic fields Hj. fdd 
is the energy of magnetostatic interactions between cells 
and £cxch describes the exchange interaction between the 
nearest neighboring cells. Boundary conditions can be 
taken into account by choosing corresponding anisotropy 
fields. 

The calculation procedure consists of the following 
steps. From (|l|) it is necessary to find the equilibrium 
direction in the yz plane for each spin (cell). This direc- 
tion is defined by the angle Oj relative to z axis. Then 
we represent the spin in its local "quantization axes" x^- , 
yj, "Zj, which are associated with the spin equilibrium 
direction (Sj || Zj): 



cosOi + S,' sine,- 
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Now we analyze small deviations for each spin from equi- 
librium. It is convenient to describe small magnetization 
oscillations by conventional spin-wave technique. Us- 
ing a linearized Holstein-Primakoff transformation [Q, 
we introduce complex variables a*, a, which are classi- 
cal analogs of creation and annihilation operators: 



S"' =-S + a*aj, ~ (oj + a*)V2S/2, 



(5) 



[aj - a*)V2S/2i 



The quadratic part of the energy (|l|) can be written in 
the form: 



1 



^ J2 ( + ^4«r + ) , (6) 



k,l 



where Afc; = A^f. and B^z = -B/fc- 

To eliminate the nondiagonal terms from we shall 
use the linear canonical transformation : 



AT 



(7) 



k=l 



where Cfc and c^. describe complex variables of normal 
spin-wave modes. From gyromagnetic equations we have 
the algebraic equations for eigenvalues (standing spin 
wave frequencies) and eigenvectors uik-, vik- 



N 

i^kUjk = ^ {Ajiuik + Bjivik), 
1=1 

N 

-uJkVjk = ^ {AjiVik + B*iuik). 
1=1 



(8) 



This system must be supplemented by the orthogonality 
and normalization conditions: 



N 



N 



~ ^ll^jk) = Slk, i.'^jlVjk-UjkVjl) = 0, 

N N 

Y^^i^^*ik ~ ""jkVik) = 6ji, ^ [ujkVik-UjkVik) = 0. (9) 



k=l 



fc=l 



The resulting quadratic energy form describes N inde- 
pendent normal modes (harmonic oscillators): 



N 



(10) 



k=l 



The evolution of each normal spin-wave mode is char- 
acterized by frequency {ujk) and relaxation rate {rjk)'- 
Ck oc exp(— iwfct — rjkt). From equipartition, white ther- 
mal noise {ksT) is applied to each independent mode. 
Using the solution for a damped harmonic oscillator with 
additive noise (see, 1^, Q), we calculate the spectral den- 
sity of the magnetic fluctuations: 



2{T]k/uJk) k^T 



^ki 



(11) 



where T is the temperature, (...) means thermal aver- 
aging. Note that the total system has thermal energy 
NkBT. 



III. APPLICATION TO A MR HEAD 

As a general application to a GMR head we take the 
z direction to be the cross track direction along which 
the current is applied and thus the deviations of mag- 
netization in the y direction are responsible for magne- 
toresistance. The power spectral density due to random 
magnetization fluctuations can be written as , (| : 

Svv{(^) = CySrnymy{(^), (12) 



Cv = hiUdR/dH,^t)idH,^t/dmy), (13) 

where 



Sm^mAuj)= I {Smy{t)5my{Q))e"^*dt. (14) 



defines the magnetization noise. Here Sniy denotes the 
net random deviation of the normalized total magnetic 
moment: 



AT 



dm,, — — 



(15) 



Using Eqs. (|^) and (^, we rewrite Smy(t) in terms 
of the normal mode variables Cfe(t) and Ci(t). Thus, the 
correlation function (n4) can be rewritten in terms of cor- 
relation functions (cfc(F)c;*(0)) and (c^(t)cz(O)). Utilizing 
(|ll|), we finally obtain: 



AT 



k=l 



Fkito) 



Wk 



[ujk - + % i^k + 1^)^ + % 



N 



(17) 



Thus the magnetization noise power ( jig ) is proportional 
to the temperature and inversely proportional to the 
film volume, and is the sum of mode spectral functions 
{Vk/<^k)Fk{(jj) with some "weights" Wk (|l7|). In the case 
of iV = 1 we obtain the result of , ||] . 

The summation of amplitudes of spin waves Ujk + Vjk 
in ( prj) qualitatively means an integration of oscillations 
in the standing mode (see, Fig.l). For example, if equi- 
librium is along the z direction {9j — 0) the integration 
over the modes in Fig. lb and Id vanish and only modes 
from Fig. la and Ic will contribute. Due to partial can- 
cellation in the odd spin- wave standing mode (see, e.g., 
Fig. Ic) the weight rapidly decreases with increasing mode 
wave number. 
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IV. EXAMPLE: TWO MODE MODEL 

Let us consider now the case of two-cell discretization 
(Fig. 2). We take into account the "easy" anisotropy field 
Hk along the z axis and the "hard" axis field AttMs per- 
pendicular to the plane. In this example only magneto- 
static interaction between cells is included taken in the 
form: (?i7)'(CSi • S2 - CySfS^-C.SfSD/V.^ where C, 
and (^z are the factors, which depend on the magne- 
tostatic energy for a given film geometry. For example, 
for a film of 4 x 180 x 250 nm^ C « 0.23, Cy « 0.17 and 
Cz ~ 0.063. 

The energy quadratic form in terms of spin deviations 
(^ can be written as: 

S^^'' /h = A{alai + a2a2) -f- i3(aja2 -I- 0201) -I- C{aia2 

+ aja2) + {D/2){aiai + a^a* -I- 0202 -I- ajOg), (18) 



where 



1 



A = 7{iJK(cos^ - - sin^ e) - Ms[C - 3(Cy sin^ 

cos^ 0)] + 2ttMs + Ho cos{0H - 0)}, 
B = jMsC + C, ~j{2nMs + Hk sin^ 0/2), 

C = -3-/Ms{Cy cos^ + Cz sin^ 0)/2. 

From (^) and ^ wc obtain: 

"11 = U21 uUi + vVi, vii = V21 = uVi + vUi, 

U12 = -U22 = UU2 + VV2, V12 = -1^22 = UV2 + VU2, 

where 



(19) 



Uk = \ : , V 



4uJk 



D A-Vt 



Vk = {-If 



k in-{-ifB-LUk _ lA + n 



2n ' 



UJk 



= J[Q-i-lfB]^ -C^, n = ^A^-\D\^. (20) 



For equilibrium in the cross track direction the two modes 
correspond to coherent and fanning about the z di- 
rection (similar to the modes illustrated in Fig. la and 
lb, respectively). From Eq.(|l9|) we can see that the 
"weight" of the second mode is always equal to zero: 
W2 oc {ui2 + V12 + U22 + V22) COS0 = 0. Thus, the two- 
mode model is reduced with some modification to the 
case of the single-mode model. 
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Figure captions. 
Fig. L 

Standing spin-wave modes (qualitatively) with odd (a, 
c) and even (b, d) numbers of oscillation humps. 
Fig.2 

Film with two-cell discretization. 
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